Resonant hyper-Raman scattering in spherical quantum dots 



00 
On 
On 

> 
O 

in 



E. Menendez-Proupin and C. Trallero-Giner 

Departamento de Fisica Teorica, Universidad de La Habana, Vedado 
10400, La Habana, Cuba 



A. Garcfa-Cristobal 
Departamento de Fisica Aplicada, Universidad 
E-4-6100 Burjassot, Spain 
(February 1, 2008) 



Valencia 



A theoretical model of resonant hyper-Raman scattering by an ensemble of spherical semiconduc- 
tor quantum dots has been developed. The electronic intermediate states are described as Wannier- 
Mott excitons in the framework of the envelope function approximation. The optical polar vibra- 
tional modes of the nanocrystallites (vibrons) and their interaction with the electronic system are 
analized with the help of a continuum model satisfying both the mechanical and electrostatic match- 
ing conditions at the interface. An explicit expression for the hyper-Raman scattering efficiency is 
derived, which is valid for incident two-photon energy close to the exciton resonances. The dipole 
selection rules for optical transitions and Fr6hlich-like exciton-lattice interaction are derived: It is 
shown that only exciton states with total angular momentum L = 0, 1 and vibrational modes with 
angular momentum l p = 1 contribute to the hyper-Raman scattering process. The associated exciton 
energies, wavefunctions, and vibron frequencies have been obtained for spherical CdSe zincblende- 
type nanocrystals, and the corresponding hyper-Raman scattering spectrum and resonance profile 
are calculated. Their dependence on the dot radius and the influence of the size distribution on 
them are also discussed. 
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I. INTRODUCTION 



The research on semiconductor quantum dots (QD's) 
has undergone a dramatic increase in recent years, stim- 
ulated by their foreseen applications in optics and elec- 
tronics technology and also due to their nonlinear opti- 
cal properties. [Q Quantum dot systems based on III- 
V materials as well as nanocrystallites of II- VI com- 
pounds embedded in glass have been thoroughly inves- 
tigated (for a review see Ref. Q ) . Among the scattering 
mechanisms present in polar semiconductors, the opti- 
cal phonon emission is known to play a dominant role 
in QD's, which can be experimentally investigated by 
employing a number of methods, such as infrared ab- 
sorption and Raman scattering. The successful in- 
terpretation of light scattering by optical phonons relies 
upon a good knowledge of the normal vibrational modes. 
In recent years a phenomenological continuum theory of 
optical phonons in nanostructures has been elaborated 
[po|— |l2|] , which is in good agreement with ab initio calcu- 
lations and allows to explain the resonant Raman scat- 
tering intensities of phonon modes induced by interface 
roughness in quantum wells. The theory has also 
been generalized to deal with quantum wires and quan- 
tum dots, and used to study resonant Raman scattering 
in these systems. |H]-|l6f In particular, the formalism is 
applicable to II- VI semiconductor nanocrystallites em- 
bedded in glass since they can have dimensions as small 
as 13 A, in which case the mechanical confinement of op- 
tical vibrational modes has strong effects (in a QD the 



concept of phonon as an excitation in a periodic system 
labeled by a wave vector is lost, and therefore we will use 
the term vibron to denote the QD vibrational modes). As 
it has been theoretically demonstrated, optical vibrons 
in spherical nanocrystallites have a mixed longitudinal- 
transverse character and can exhibit different symme- 
tries. |uj Raman scattering provides a useful tool to in- 
vestigate experimentally these vibrons, but only spheri- 
cally symmetric modes are accessible to this technique. 
In the search of complementary experimental techniques 
which can overcome this limitation, hyper-Raman (HR) 
spectroscopy appears to be a suitable candidate. |l~7|] Re- 
cently, HR spectroscopy has been used to study the opti- 
cal vibrational modes of CdS and CuBr QD's. |l8|-]20|l In 
this paper we present a theoretical model which allows to 
study hyper-Raman scattering (HRS) by optical vibrons 
under resonance conditions, and illustrate it by perform- 
ing numerical calculations of the scattered intensities in 
CdSe nanocrystallites. As far as the exciton-lattice in- 
teraction is concerned only the Frohlich coupling is con- 
sidered here, though it is admitted that the deformation 
potential interaction might be of importance for photon 
energies far from the excitonic resonances. The selection 
rules for the scattering process are worked out and it is 
shown that in fact HRS can be used to probe optical 
vibrations with non-spherical symmetry. 

The paper is organized as follows. The main concepts 
underlying the hyper-Raman scattering and the descrip- 
tion of the nanocrystal vibrational modes and excitons 
are outlined in Sec. |n|. Section III contains the theoret- 
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ical expressions for the matrix elements and HRS effi- 
ciency, which are used to analyze the selection rules of 
the scattering process. The calculation of the exciton 
and vibrational spectra of CdSe nanocrystallites are pre- 
sented in Sec. [j^, and the numerical results obtained for 
the HRS efficiency are reported and discussed in Sec. [v|. 
Section [v| is devoted to the main conclusions of the work 
and final comments. 



where V is the normalization volume, r/i (i — 1,2) and 
rj s are the nanocrystal refraction indices at w; and to s , 
respectively, and c is the velocity of light in vacuum. Ac- 
cording to the Fermi's Golden Rule, W can be written 
as 



W((Ji,Wj,tJ B ;R) = 



M p (cji,LOj,Lo s ;R)\" 



xS(uji + ujj — uj p — <jj s ) 



(3) 



II. THEORY 

A. Hyper-Raman efficiency 

Hyper-Raman scattering is a nonlinear process which 
consists of the absorption of two photons of frequency 
u>i, wave vector k^, and polarization (i = 1,2), and 
the emission of one photon (w s ,k s ,e s ) with the simulta- 
neous excitation of a number of vibrational modes. H ] 
Due to its nonlinear nature, it is convenient to express 
the HR yield by the normalized (intensity-independent) 
scattering cross section [E2| 
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1 d 2 Pjj 
Iilj duj s dQ s 



d Pi- 



(1) 



the 



Ii and Ij being the excitation intensities and duj df2 
scattered power per unit of solid angle fl s and unit fre- 
quency. The scattered power can be related to the proba- 
bility per unit time W(u>i, ujj,Lu s ; R) of a single scattering 
event in a QD of radius R, so that we can express the 
normalized differential cross section as [E2l 
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W(uJi,u!j,u! s ;R) 
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where M p is the scattering amplitude. Each scattering 
event is accompanied by the emission of a vibron of fre- 
quency uj p , and the total scattering probability is ob- 
tained after summing over all possible vibron states p. 
Only vibron emission (Stokes) processes will be consid- 
ered in this work. 

The delta function in (||) can be eventually replaced 
by a Lorentzian, 



A(uj s ) = - 



5(u>i + bJj —Up~ UJ S ) 



7T (flU>i + htUj — TlLiJp 



(4) 



in order to take into account the effect of the vibron life- 
time Tp = h/Tp on the line broadening of the HR spec- 
trum. 

The scattering amplitude in (|^) can be calculated in 
fourth-order perturbation theory. Under resonance con- 
ditions the important contributions to M p are illustrated 
by the diagrams in Fig. 1. In Fig. 1(a), the exciton cre- 
ated in the state after the absorption of a photon 
(a;,-,kj,ej) is first scattered to the state I//2) by the ab- 
sorption of a second photon (a;,-, kj, e^). In the next step, 
the interaction with the lattice induces an excitonic tran- 
sition from I/Z2) to I/X3), accompanied by the creation of 
a vibron of frequency uj p . The exciton finally recom- 
bines emitting a photon (w s , k s , e s ). The amplitude cor- 
responding to this process is 



M (a) 



(F\H^_ R \^ 3 ) (mIHe-M (»i\H%L r \I) 



y - 



E„ 



(5) 



where \I) (\F)) is the initial (final) state of the scat- 
tering process, E /li and T /li are the energies and lifetime 
broadenings of the excited electronic states |//j) in the 
QD. He-l and He-r are the Hamiltonian operators for 
the interaction of the electronic system with the lattice 
and radiation field, respectively. 

Since the interesting range of hui+hiVj for resonant HR 
spectroscopy occurs around the fundamental absorption 



edge Eo, the relevant resonances in the matrix element 
(0) will occur at energies 

huji + Tiijj — Ep 2 (incoming resonance) , (6a) 
hu> s = Ep 3 (outgoing resonance) . (6b) 

For the process shown in Fig. 1(b) the scattering ampli- 
tude is 

(7) 



(F|g^ fl | M3 )( M 3|g^|^)(^2|gi ; -L|/Xl)(/X 1 |g^ fi |/) 

(huj s - Ep 3 + iTp 3 )(huji - hujp - Ep 2 + iTp 2 )(fiuji - E^ + iT^) 



rU) 



(i) 
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In Eqs. (jsj) and ([?]) it is implicitly understood that 
the sum over non-topologically equivalent diagrams (i<->j) 
must be taken into account. The examination of the 
energy denominators in Eq. (0) clearly indicates that 

\Mp | -c \Mp a ^ | in the resonance region huji + Tiujj ~ Eq. 

Therefore the contribution from M p has been dropped 
out from our theoretical model. 

Finally, it must be noted that in real samples the 
nanocrystallites present a distribution over size and 
shape. We intend to study here HRS by an ensem- 
ble of spherical QD's characterized by a distribution 
over radii F(R). The corresponding average (intensity- 
independent) HRS efficiency [p3| is given by 

\ / ~ (V D ) J h dw s dn s F ( R > dR ' « M 

where (Vb) is the average quantum dot volume. 

B. Exciton states 

The intermediate electronic virtual states appearing in 
the HR process ( see Eq. (^) ) are taken to be size-confined 
Wannier-Mott excitons, treated in the framework of the 
envelope function approximation (EFA). The details of 
this approach can be found in Ref. JTr| , where it is em- 
ployed for the study of resonant Raman scattering in 
QD's. We summarize here, for notation purposes, the 
main expressions of the exciton model. In the absence of 
electron-hole interaction, the wave function of an uncor- 
related electron- hole pair (EHP) in a spherical QD can 
be written as 

$a,£,Af(r e ,rh) = ^2 (lelhm e m h \LM) 

(f>n e ,l B ,m e (re) <f>n h ,l h ,m h (rh) , (9) 

where a is an abbreviated notation for a = (n e , nu, le,lh)- 
The single-particle (electron or hole) wave functions are 
given by 

(r) = RnAr)Yi, m {^) ■ (10) 

Yi,m{^) (I = 0,1,... and m = —I,..., I) are the spheri- 
cal harmonics J24|, and R n ,i(r) is the radial part of the 
wave function, which depends on the confinement poten- 
tial. The corresponding energy is denoted by E n j. L 
and M are the quantum numbers corresponding to the 
squared total angular momentum, L 2 , and its projection 
along the z-axis, L z , respectively, and (L lh m e rrih\L M) 
are the Clebsch-Gordan coefficients. [P5[ 

The eigenfunctions of the full Hamiltonian includ- 
ing the electron-hole Coulomb interaction (exciton wave 
functions) are expanded in terms of the EHP states (|^), 

* M (r e ,rfe) = &N,L,M,p(r e ,Th) = 

^ Cn,l,m,p(u) ®a,L,M(r e , r h ) , (11) 



where N and P are the quantum numbers correspond- 
ing to the energy and inversion operator, respectively. 
Note that the state defined by Eq. ( |TT| ) has definite par- 
ity: P = 1 (even parity) if l e + lh is even and P = — 1 
(odd parity) if l e + lh is odd. By applying the matrix di- 
agonalization technique described in Ref. jig] we obtain 
the exciton energies as the eigenvalues of the secular 
equation and the coefficients C M (a) as the corresponding 
eigenvectors. 



C. Exciton-lattice interaction 

Excitons in semiconductors interact with lattice vibra- 
tions by means of deformation potential and piezoelec- 
tric interactions. In polar materials, the long-range elec- 
trostatic field associated to the optical vibrations intro- 
duces a different coupling mechanism, the Frohlich inter- 
action. It is now well established that, despite its dipole- 
forbidden character |23]j , the Frohlich interaction plays 
an important role in one-phonon Raman scattering by 
bulk zincblende semiconductors. |^6| On the other hand, 
in systems which lack translational invariance, like QD's, 
Frohlich-induced Raman scattering becomes allowed and 
therefore it is expected to have a dominant role. Ac- 
cordingly, in our model for HRS we consider only the 
Frohlich- like interaction between excitons and vibrons. 

In order to describe the polar optical vibrational modes 
(vibrons) of spherical QD's we rely upon the results of 
Ref. E3| , where a macroscopic continuum model coupling 
the mechanical displacement and the electrostatic poten- 
tial is developed. For the case of complete confinement it 
has been shown that the normal modes with non vanish- 
ing electrostatic potential exhibit a mixed longitudinal- 
transverse character. They are labeled by a set of integer 
numbers p = (n p , l p , m p ), which are related to their sym- 
metry properties. More details on the calculation of the 
displacement, electrostatic potential and frequency u>n r .i r 
associated to these modes can be found in Refs. |1J] and 
II- 

The exciton-vibron interaction Hamiltonian operator 
can be written as 

He-l = e<£_F(r e ) - e<p F (r h ) , (12) 
where — e (e > 0) is the electron charge and 



eip F (r) = 




x(Y lptmp (il)b npJpimp +H.c) . (13) 

Here b n i m is the vibron annihilation operator and 
H.c. means Hermitian conjugate. The explicit form of 
the radial function &n p ,l p ( r ) can 

be found in Refs. @,||. 
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Since uj np ,i p ~ (wl is the bulk LO-phonon fre- 
quency at the T-point), we will omit hereafter the factor 
^/w £ /w n j in (p^|). The Frohlich constant Cp is given 
by 



l2Tre 2 hcj L ( — - 



£o0 



(14) 



e aoo and e a o being the high- and low-frequency dielectric 
constants of the crystallite material. 



III. MATRIX ELEMENTS AND SELECTION 
RULES 

In this section we analyze in detail the matrix elements 
appearing in Eq. (|5|) , and derive from them the selection 
rules for the exciton and vibron states which participate 
in the HR process. 

The matrix elements (fii\Hg_ R \I) and (F\H^ R \/j, 3 )) 
for direct allowed optical transitions between valence (v) 
and conduction (c) bands are given by J16| 



fai\H$_ R \I) 




( e i ' Pcv) ffii 



(15) 



where p cv is the interband momentum matrix element 
between valence and conduction Bloch functions at k = 
0, and mo is the free-electron mass. A similar expression 
can be obtained for (^\H^_ R \F). The exciton overlap 
integral / M is given by 



Tn,l,m.p 



N.L.M.P 



>,r) d 3 r = S L nS 



M,o 



xS P ,i X^ 1 )' + 1 c N,o,o,i(n e ,n h , I, I) 



n e ,n h I 



R n ^i{r)R nh ^{r)r 2 dr . (16) 



Hence, the annihilation of the first incoming photon 
with frequency uii creates an exciton in the state with 
zero angular momentum and even parity, 

/i 1 = (JV 1 ,£ 1 =0,Mi = J Pi=l) , 



and analogously, the scattered photon of frequency lo s 
is emitted upon the recombination of an exciton in the 
state 

P3 = (JV 3> L3 = 0,M3 = 0,P3 = 1) ■ 

Both exciton states have L = and P = 1 because the 
interband transitions induced by the radiation field re- 
quire the electron and hole to have equal orbital angular 
momentum quantum numbers, l e — In- 

A different situation is found when considering the 
scattering between exciton states induced by the sec- 
ond incoming photon (matrix element (/i2|i?^l B |/ii) in 
Eq. (j^)). Two possibilities arise in this case p7[ : On 
the one hand there may be transitions accompanied by 
scattering of the electron (or hole) between QD levels cor- 
responding to different bands (interband excitonic transi- 
tions) , and by the other hand there are also transitions in 
which the electron and hole remain in the same band, and 
only the exciton envelope function is affected by the inter- 
action with the radiation field (intraband excitonic tran- 
sitions). In III-V and II- VI compounds, and for Tiu>i+TuOj 
around the fundamental absorption edge, the interband 
excitonic transitions are due to remote bands which are 
far enough in energy, thus giving a negligible contribution 

to the scattering amplitude Mp . Therefore we consider 
only intraband excitonic transitions, whose matrix ele- 
ment can be written, in the dipole approximation, |27j] 
as 



(fi2\H U) 



K 2 (r e ,r h ) 



H e-R ~ H h-R J *w r h) d i v e d A r h 



r(J) 



where 



n v-R — 



1 



2tt77, 



m 



(17) 



(18) 



p, y = — ih Vj/ is the linear momentum operator, and m e 
(m/j) is the electron (hole) effective mass (taken to be 
positive) . 

Using the expansion (|ll]) in ( pj] ) we get 



(m\H ( g_ R \ni)=e- 



C N 2 ,L 2 ,M 2 ,P 2 ( a ') C NxM,Mx,Pii a ) { a',L 2 ,M 2 



V J J a,a f 



e,- • Pe 



Ph 



a, Li, Mi 



(19) 



Let us now concentrate on the electron part of the matrix element appearing in (|19|). By making use of the operator 



identity p e — (im e /H) 



H e ,r e 



where H e is the single-electron Hamiltonian, it is easy to find that 
e • p e 



n' e ,n' h ,l' e ,l' h ,L2,M-2 



n e ,n h ,l e ,l h ,Li,Mi ) =5, 



E nc j e )(n' e ,n' h , l' e , l' h , L 2 ,M 2 \e ■ r e \n e , n h , l e , l h , Li, Mi) 



(20) 
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The matrix element in the previous equation can be pendix. Moreover, for L\ = Mi = (l e = lh), Eq. Q2C 
greatly simplified by taking advantage of the theory of reduces to: 
angular momentum, as it is explicitly shown in the Ap- 



n e ,n h ,l e ,l h ,L 2 ,M 2 



iR 1 , 



e • p £ 



n e , n h ,l e , l e , 0, ) — S n ' htflh Si> h ,i h <>l 2 ,i Sp 2 ,-i a >M 2 



le + l 



2l e + 1 



Sr 



l'.,le + l 



I, 



2L + 1 



(21) 



r -1 r 

R Jo 



Rn',i'(r) R n ,i(r) r 3 dr 



(22) 



The final expression for the matrix element (|19|) can be arranged in the form 



\M2 We-R Ml/ — T~S 7f7\ 2 

TiR yV W Ujifj 



f( e ) _ f{ h ) 



(23) 



where we have introduced the energy Eq — Ti 2 /2moR 2 {R is the average QD radius) to make T^ e ^ dimensionless. 
The explicit expression for T^}^^ is 



RY (E n 



R 



^jfl^lM2 = S L2-1 $P 2 -1 a M 2 5 iaJl' h ,lh S n' h ,n h C* 2 {a') G ^ (a) 

a,a f 

f>t' Z e + 1 — \ / „, C , -i Sv i e -i G ne j e ^> n > i> 



E, 



le + l 
2/ e + 1 



2l e + 1 



(24) 



A similar expression holds for JF^L M2 after the ex- 
change of the subscripts e and h. 

An important consequence to be drawn from (|24|) is 
that after the absorption of the second incoming photon, 
the exciton must be in the state: 

At2 = (jV 2 ,L 2 = l,Af 2 = 0,±l,P 2 = -l) . 

Otherwise stated, two-photon absorption generate exci- 
tons in L = 1 states, in contrast to one-photon tran- 
sitions for which the final exciton state is necessarily 
L = 0. Moreover, the factor <5/',;±i in Eq. ( p4| ) intro- 
duce the parity selection rule Pi = 1 — * P 2 = — 1, 
indicating that the intraband transitions are accompa- 

C F 



nied by a change in the parity of the excitonic state. If 
the incident light is linearly polarized with ej || z then 
a M 2 = <5m 2 ,o (see the Appendix), whereas for circular po- 
larization ej || (x ± iy)/y/2 we have om 2 = <5a/ 2 .±i ( x, y, 
and z represent some system of orthogonal axes attached 
to the laboratory frame). Hence, the L 2 = 1 excitonic 
states participating in the HRS process have M 2 = 
for linearly polarized light and M 2 = ±1 for circularly 
polarized light. 

Now we turn to the matrix element of the exciton- 
lattice interaction UjAHe-l\^2}' If we select the vibron 
creation terms in (|l^) and make use of (pi]), we get the 
expression 



(^■a\He-l\^2) - — ^ ^2 G* N3 x 3 M 3 ,P 3 { a ")GN2,L 2 ,M2,P2{ a ') 
a' ,a" 

]T (a",L 3 ,M 3 $„ piip (r e )F^ mp (r! e )-$ np , p (^)y^ mp (^) a',L 2 ,M 2 ) 



(25) 



which consists again of separated electron and hole 
contributions, each of them being a sum of matrix el- 
ements of spherical tensors. Therefore the method used 
previously for the intraband exciton-photon matrix ele- 
ments (see the Appendix) applies also here. By following 
that procedure, and considering that L% = 1, L3 = M 3 = 



0, Eq. ( |25| ) is reduced to 
(/i 3 |-ffs-i|M2) = ^= 
where 



-rj(e) _ <yj(h) 



(26) 
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n'A'-^n'JA'' ' 



(27) 



n' ,V — >n" ,V' 



Rn",i" (r) Rn>,i> (r) ®n p ,i p (r) r 2 dr 



(28) 



Ah) 



after 



An analogous expression to (|27j) holds for H^- 
the exchange of the subscripts e and h. 

It follows from (^) that the quantum numbers of the 
emitted vibron are fixed to be l p = 1 and m p = M%- 
The different states m p can be discriminated by select- 
ing adequately the polarization of the exciting light: 



M 2 = 



= if e, 



and M 2 = 



±1 if 



Gj || (x±iy)/y/2. Thus, the HR selection rules in a spher- 
ical QD can be expressed schematically by the following 
sequence of states. 



(Li = 0, Mi = 0) — ► (L a 

(Z p = l,m p =M 2 ) 



l,M a = 0,±1) 
(L 3 = 0,M 3 = 0) , (29) 



which contrasts with that corresponding to Raman scat- 
tering HP] 

(L 1 =0,M 1 = 0) ( ' P= °^ P=0) (L 2 = 0,M 2 = 0) , (30) 

where only l p = vibrons can be excited. The compar- 
ison between ( p9| ) and ([30j) makes up clear the complc- 
mentariness of both light scattering techniques to study 



the vibrational spectra of spherical QD's. In particular, 
we have shown that hyper-Raman spectroscopy can be 
used to investigate the l p = 1 vibronic states, not ob- 
servable in Raman scattering. In addition, the analysis 
of the corresponding resonance profiles (scattering effi- 
ciency vs. hu>i + fuvj) can be useful to reveal L = 1 ex- 
citon states. Of course, the selection rules discussed can 
be relaxed when going beyond our simplified treatment 
of the HR process, e.g. including valence band mixing. 
Nevertheless, for R smaller than the bulk exciton Bohr 
radius the separation in energy between the hh and Ih 
levels induced by the confinement should lead to a small 
amount of hh-lh admixture. 

Finally, by inserting Eqs. (g§), @, and @ into (|), 
we obtain the following compact expression for the nor- 
malized scattering cross section 



1 d 2 a ij 



R 



Ii dwsdfls \RJ <q irj . y UJj 

|2 



Vs 



\M np (u)i,Wj,w s ;R)\ A(uj s ) , 
where M. n is the dimensionless amplitude 



(31) 



Mr 



En 



E 



f * 

J «3 



o/(e) 



(<0 



M2 Ml- 



(hu s - E P3 + iT P3 )(huji + Tiujj - E^ + iT^fiuji - E^ + il\ 



(32) 



and 



(TO = 27T- 



R 2 



|e» ■ p c 



|e s • p c 



r 2 



TTin C 5 



m 2 E 2 



E 2 R 



(33) 



The expression (|31| ) is suitable for the calculation of 
the scattered spectrum as a function of the HR shift, 
lu s — cui — ujj . If we are interested in the resonance behav- 
ior of the scattered intensity we can obtain the average 
HRS efficiency by integrating Eq. (B) over ui s , 



Eq (~ 10 meV), and the parameters of Table | a value of 
10~ 4 cmMW^sr" 1 is estimated for (<r )/ (V D ). In the 
case of an ensemble of QD's the HRS efficiency (|34|) is 
replaced by the average over the QD orientations with 
the angle arccos^ • e s ) fixed. 

From now on we focus on the degenerate case, cj, = 
u>j = u>i, with 2loi in the region around the excitonic 
transitions, which is the usual situation for resonant HRS 
experiments. 



dfl s 



(Vi 



D, 



E — 



Wi + wj - uj 1lp (R) 



IV. EXCITON ENERGIES AND VIBRONS IN 
CDSE NANOCRYSTALLITES 



|M„,(w i> a; J ,w i + a; J -w nj> (J2);iJ)| F(R)dR, (34) 

In the following we consider for the radial part of the 
where (<to) is the average of <7o over the QD orientations single-particle wave functions R n ,i(r) the solutions of the 
in the ensemble. For typical values of R (~ 20 A) and, infinite barrier spherical well problem 



G 




1 



*(-S?h) 



(35) 



where ji is the spherical Bessel function and the 
roots of ji(xn ) — 0. The single-particle energies are 
given by 



E n ,i — 



2m R 2 



-0)' 



(36) 



First, we present in Fig. g(a) the exciton energies, calcu- 
lated as explained in Sec. for a CdSe quantum dot as 
a function of the radius. The material parameters used 
as input are shown in Table [| Solid and dashed lines 
correspond to L — and L = 1 excitons, respectively. 

According to Ref. || we have taken into account the 
penetration of the wave function in the glass matrix by 
using an effective radius R ef = 3.367 + 0.9078i? (in A), 
where R is the nominal QD radius. 

Let us focus now the attention on the vibrational 
modes. Since the active TO-phonon branch of hexagonal 
CdSe seems to be flat the parameter (3t appearing 
in the isotropic model for the optical vibrations has been 
taken as the limit f$\ — > + , with a negative bulk disper- 
sion relation J 1 — ui\ — 0^g 2 ■ If one attempts to solve 



the vibron model of Ref. [[L4| with (3t = 0, it is not pos- 
sible to fullfil all the matching boundary conditions at 
the interface, What is expected on physical grounds 
is that the TO component of the vibron vector displace- 
ment, at forbidden frequency for pure TO phonons, is a 
vanishing or rapidly decaying function. This can be ver- 
ified in the solution for the vibrational amplitude u(r) 
(Eqs. (35) and (36) of Ref . ). The TO component of 
u vanishes for all r < R in the limit /3f, — > + . 

Figure ||(b) illustrates the allowed frequencies for the 
l p = 1 optical vibrons as a function of the QD ra- 
dius. Below lul = 213 cm -1 it can be seen the con- 
fined LO modes with both longitudinal and transverse 
components, including a surface mode contribution. |2{J 
We can observe some bending in the dispersion around 
185 cm -1 , which is identified as the Frohlich mode. ]3C| ] 
This mode arises in the framework of dielectric models, 
which neglect the effect of mechanical boundary con- 
ditions and consider flat bands for TO and LO bulk 
phonons. p9[ | The Frohlich frequency ojp is related only 
to the dielectric constants of the constituent media: lo f = 
Wy(e o +2e(, oo )/(e a0O +36600). A strong electrostatic con- 
tribution is expected in the dispersion for vibron fre- 
quencies close to ujp. In Fig. || we depict the electro- 
static potential i(r) as a function of r for QD radii 
R = 11.5 A, 16.2 A, and 21 A. At these radii, u) n \ equals 
u>f for n p — 2, 3, and 4. It can be seen that the electro- 
static potential <E>„ t \ shows an enhancement close to the 
interface whenever uj n 1 coincides with top. 



V. SCATTERING INTENSITIES 



In this Section the results obtained with Eqs. (J3l|)-(|34|) 
for the HRS spectrum and resonance profile in CdSe QD's 
are shown. The incoming (outgoing) resonances will be 
denoted by the combination of the letter I (O) and the 
exciton quantum number N of the corresponding reso- 
nant levels. It is important to remind the selection rules 
already discussed in Sec. [nj The relevant incoming res- 
onances take place at 2Tiu>i = -E/v,l=i,a/.-i whereas the 
outgoing ones appear at hui s = -E/v,l=o,o,i- 

In the next discussion we first analyze the hyper- 
Raman process for the case in which all the nanocrystal- 
lites in the sample have the same radius R. Typical HRS 
spectra are shown in Fig. ^ for different QD radii. The in- 
cident photon energy is such that fru> s — 2hLUi — hoji t i is in 
resonance with the N = 1, L = excitonic level. We have 
included in the calculation the lowest 11 QD excitonic 
levels (N = 1 to 11, for each value L = 0, 1). A lifetime 
broadening T M = 5 meV was assumed for all excitonic 
transitions. The spectra are broadened to have a full 
width at half maximum (FWHM) 2r„ p = 4 cm -1 . It is 
systematically found that the main peak of the spectrum 
corresponds to the creation of the l p = l,n p = 1 vibron, 
with a small contribution coming from the n p = 2, 3, 4 vi- 
brons. This is related to the fact that the matrix element 



n 



(e,h) 



drops off rapidly as n p increases. It must be 



noted that in the spectra 4(a) and 4(b) the n p = 3 peak 
is stronger than the n p — 2, 4 ones, while the largest con- 
tribution in spectrum 4(c) (aside from the main line) is 
due to the n p — 4 peak. This difference can be explained 
in terms of the electrostatic effects discussed above: If 
one looks back to the vibron frequencies as a function 
of the QD radius in Fig. §(b), it can be realized that 
for R = 16 A and 18 A the vibron frequency W3.1 is 
around the Frohlich frequency uip , and thus the exciton- 
vibron interaction is dominant for these modes in the 
corresponding spectra 4(a) and 4(b). For R = 21 A, on 
the other hand, tt>4,i is approximately equal to ujp and 
therefore its potential gives an enhanced contribution to 
the exciton- vibron interaction (see also Fig. ^|), which is 
reflected in the spectrum 4(c). 

The HR resonance profile (scattering efficiency vs. 
Ihwi) for the n p = \.l p = 1 vibron peak is displayed 
as a solid line in Fig. [| for QD radius R — 18 A. The 
most important feature of Fig. || is that the 01 (L = 0) 
resonance takes place at lower energy than the II (L = 1) 
resonance, which is a natural consequence from the spac- 
ing between excitonic levels being much larger than the 
vibron frequencies. It is worth pointing out that the op- 
posite situation is usually encountered in bulk semicon- 
ductors. |2^] We have also calculated the HRS intensity 
taking the Coulomb interaction equal to zero, recovering 
the free electron-hole model |l5| (dashed line), and treat- 
ing the Coulomb interaction just in first order perturba- 
tion theory (dotted line). It is apparent the exciton red- 
shift when the electron-hole interaction is included. Also, 
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when comparing the absolute values of the scattering in- 
tensities for the different approaches we see that the full 
calculation is extremely well approximated by the per- 
turbative approach and gives only slightly larger values 
than the free electron-hole model (mainly due to the en- 
hanced oscillator strengths of L = excitons). Concern- 
ing to the HR cross section the free electron-hole model 
presents identical lineshape to those displayed in Fig. ||, 
whenever 2Huji is rescaled to set the equivalent resonance 
conditions. Thus, we conclude that the excitonic effects 
on the HR resonance profile (and also on the HR spec- 
trum) of quantum dots in the strong confinement regime 
stand mostly to renormalize the resonance energies. 

Let us finally discuss the effects of the size dispersion 
of the crystallites on the HRS spectrum. We have con- 
sidered an ensemble of QD's described by a Gaussian 
distribution function F(R) centered at the mean radius 
R and with FWHM equal to 40%. For a given incom- 
ing photon energy hui, Eqs. (JsJ) determine a set of res- 
onance radii {R r } and their corresponding resonant ex- 
citon levels {E r } with lifetime broadenings {r r }. These 
are marked by arrows in Fig. |^(a), at Ihioi — 2.771 eV. 
As R departs from R r , the resonance condition ceases to 
hold rapidly. The range 5R r around the resonance ra- 
dius R r which contributes to the integral in Eq. ( |34| ) is 
typically of the order of some tens of Angstrom. In such 
a small interval, the matrix elements and the vibron fre- 
quencies in Eq. (|3^) can be considered as constants. On 
these grounds we have employed the following approxi- 
mation 

r 

An estimation for 8R r can be obtained from 

5R r = ^E7T ■ (38) 
dR \R r 

Note that as the HR spectrum evaluated at the reso- 
nance radius R r is proportional to T~ 2 and SR r is pro- 
portional to r r , the contribution of the i? r -resonant QD's 
to the average HR spectrum is proportional to 1/T r . 

In Fig. U we show the averaged HR spectrum of the 
ensemble of CdSe spherical nanocrystallites with a mean 
radius of 21 A, obtained when 2hu>i = 2.771 eV. This two- 
photon energy determines a set of incoming and outgoing 
resonances with different exciton levels for QD's with dif- 
ferent radii R r , which are listed in Table H (for the emis- 
sion of the n p = l,l p = 1 vibron). The contributions of 
the various R r are also shown in Fig. ^. The resonances 
due to other exciton levels either are too weak or are at- 
tenuated by the size distribution function. Two impor- 
tant features of the HR spectrum deserve special atten- 
tion. The first one is the main peak at about 208 cm" 1 , 
which is originated by the emission of the n p — 1, l p — 1 
vibrons in QD's with different resonance radii. The most 
important contribution is that of the outgoing resonance 



03 which takes place for R = 28.0 A. This resonance ap- 
pears so strongly because a double resonance condition 
conecting the levels N = 3. L = and N = 3, L = 1 
is almost fulfilled (see Fig. ||(a)). In addition the res- 
onance 13, associated to dots of radius R — 27.8 A, is 
very close in energy to 03. The second important fea- 
ture is the broad structure between 180 and 190 cm" 1 . 
This structure is caused by the emission of interface-like 
vibrons from all resonant QD's. Its irregular shape is a 
direct consequence of the vibron frequency dispersion as 
a function of QD radius. 

Figure |?j shows the variation of the HR lineshape with 
small variations of the incoming photon energy for an 
ensemble of QD's with a mean radius of 18 A. When 
2fiuoi — 2.771 meV increases (decreases) by 100 meV, the 
radii of resonance get smaller (larger) by about 1 A (see 
Fig. ||(a)). Hence, the principal peak of the spectrum is 
shifted to high (low) frequencies, as can be seen in the 
four lower curves of the figure. Nevertheless, an oppo- 
site effect appears in the upper spectrum (change from 
2hwi = 2.771 eV to 2%uj 1 = 2.971 eV): In this special sit- 
uation the higher resonance radii get closer to the mean 
radius R — 18 A and become dominant in the spectrum. 
This effect is equivalent to an increase of the mean radius. 
It is important to note that the FWHM of the spectrum 
at 2hu>i = 2.771 eV is larger than at 2fiwi = 2.671 eV as a 
result of the competition between the large contributions 
of the resonance radii. 

Finally, the dependence of the HR lineshape on the 
mean radius for 2hwi = 2.771 eV is shown in Fig. H. The 
main features can be explained with similar arguments 
as in the previous figures. For R = 15 and 16 A the 
lineshapes are practically identical due to the resonance 
01 in the dots of radius R = 18 A. For R = 18 A the 
maximum is slightly shifted to high frequencies reflecting 
the influence of other resonances associated with larger 
radii. For R = 21 A the main peak is not Lorentzian, 
showing the contribution of several equally-strong reso- 
nances. And for R > 21 A the resonance 03 becomes 
dominant. 



VI. CONCLUSIONS 

We have developed a formalism to calculate the nor- 
malized hyper-Raman scattering efficiency in spheri- 
cal semiconductor quantum dots, considering confined 
Wannier-Mott excitons as the intermediate states. The 
exciton-lattice interaction is assumed to occur via 
Frohlich-type coupling. It has been shown that hyper- 
Raman spectroscopy can be used to probe the l p = 1 
vibrational modes. In addition, each particular mode 
(m p = 0,±1) can be selected by properly choosing the 
polarization of the incident light. With linearly polar- 
ized incident light, only m p = vibrons contribute to 
the scattering whereas by employing circularly polarized 
light the m p — ±1 modes are active. The details of 
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the calculations of the polar optical vibrational modes in 
CdSe QD's have been discussed and the eigenfrequencies 
and electrostatic potentials of the l p = 1 modes have been 
presented as a function of the nanocrystallite radius. It 
has been demonstrated that when their frequencies are 
close to the Frohlich frequency, they undergo an increase 
in their electrostatic surface character. 

The calculations of the hyper-Raman spectra show 
that the most prominent peak is due to the emission of 
n p = l,l p = 1 vibrons. The other contributions are due 
to the interface-like l p = 1 vibrons. The presence of 
surface electronic excitations could increase their role in 
hyper-Raman scattering. 

The effect of the electron-hole Coulomb interaction has 
been found to be of limited importance for the hyper- 
Raman scattering in the strong confinement limit rele- 
vant to the QD's analyzed here. The absolute values of 
the HRS intensity are slightly enhanced by the Coulomb 
interaction as compared to the free electron-hole model 
values. The main excitonic effect seems to be the renor- 
malization of the resonance energies. 

The size dispersion of the nanocrystallites is shown to 
give rise to a complex behavior of the main peak position 
in the HR spectrum as a function of the laser energy and 
the details of the distribution over radii. An additional 
effect is the distribution of the signal due to interface-like 
vibrons over a broad band of frequencies from 180 to 190 
cm -1 . 
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APPENDIX: REDUCTION OF THE MATRIX 
ELEMENTS 

In this Appendix we show in detail how to calculate 
matrix elements of the form 



(n' e ,n' h ,l' e J' h ,L',M'\T kq \n e ,n h ,l e J h: L,M) , (Al) 

where T kq is the qtYi component of a fcth-order spherical 
tensor (single-particle) operator T k , \n e ,nh,l e ,lh, L, M) 
represent an electron- hole pair state ( see Eq. (^) ) , and 
(_L,M) are the corresponding total angular momentum 
quantum numbers. 

First of all, we factorize the dependence on M' , M and 
q by applying the Wigner-Eckart theorem: 

«, n' h , l' e , l' hl L', M'\T kq \n e , n h , l e , l h , L, M) 



= (-1) 



L'-M' 



V k L 
-M' q M 



V2U + 1 



(n' e , n' h , l' e , l' h , L'\\T k \\n e , n h ,l e , l h , L) 



(A2) 



The factor (. . . , L'\\T k \\ . . . , L) does not depend on q and 
is called reduced matrix element. It can be further sim- 
plified by taking into account that T kq is a single-particle 
operator. Let us suppose, without loss of generality, that 
this operator only acts on the electron coordinates. Then 
the following reduction formula can be applied [f25[ 



(n' e , n' h , l' e , l' h , L'\\Tj, e) \\n e , n h , l e ,l h , L) 
X , fi,,, (-\) k + l >>+ L + l e J L ' L k 

*n' h ,n h \l h { J-J | le j/ ^ 

x V(2i+1)(2^ + 1)«, l' e \\T^\\n e , l e ) 



(A3) 



Finally, by inserting (A3) into (A2) we ac hieve the com- 
plete simplification of the matrix element (Al). 

Let us illustrate this procedure taking as example the 
matrix element in (^0|) . As a previous step we must real- 
ize that the operator e • r e can be expressed in terms of 
spherical irreducible tensors, 



<?=-! 



with the definitions 



1'/ 



— r e Y ljq (Q e ), (A4) 



a 



a±i = 



V2 



(A5) 



Now, by applying Eqs. ( |A2j ) and ( [A3j ) we obtain 

«, n' h , l' e ,l' h , L 2 , M 2 |e • r e |n e , n h , l e , l h , L\,M\) = S <>nh (-1)1+^+^1+^-^ 

(La 1 L x \ / L 2 L x 1 \ / l' e 

M2 ~ Ml ^ -M a M 2 - Mi Mi )\l e l' e l h j\0 

x V(2L 2 + l)(2Li + l)(2l' e + l)(2/ e + 1) (^J R K ,vSr) Rn.,l a (r) r 3 dr" 



(6) 



Proceeding in a similar way it can be shown that an used the following reduced matrix elements 
expression analogous to (^) (with subscripts e and h ex- 
changed everywhere) holds for the matrix element of the 
hole operator e-rv To obtain the last expression we have 
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; / ,inm ii i\ ( f r, f \ r. / \ 3 , \ ,,/inrim Moreover, if we particularize (til) to the case Li = Mi = 

(nJWT^nJ) = ^y{J ^'A^RnMr dr j {I \\Yi\\l) „ ^ = ^ and eyaluate the 3^ Qj symbolg we get 

the result shown in Eq. (pl|). 

Applying the same method the following expression is 
ni\w iu\ / l)(2fc+ X) f I' k I \ /cA obtained for the matrix elements of the electron-lattice 

(JIIYfeHZ) = (-1) y ^ Q Q Q j . (8) interaction 

C F 



( n yilXL 3 M 3 \eMre)\ny h l'J' h L2M2) = ^.^^.^-^^^^^^^-^^(Z^ + 1)(2L 3 + 1) 

£ \M + w + m + ( i S ) ( I o p f ) { § ff 5[ } ■ ( 9 > 



In this expression it is implicit that only the cre- 
ation part of ipp is acting on the ket. Replacing L2 = 
1,L 3 = M 3 = 0, the Eq. (g) is reduced to Eq. @ and 



n e >^e — m e'^e 



is defined in Eq. 
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FIG. 1. Feynman diagrams which give the main contribution to the resonant hyper- Raman scattering. 



FIG. 2. (a) Energy levels of N = 1, ...7 and L — (solid lines) and L = 1 (dashed lines) excitons in a spherical CdSe 
nanocrystallite as a function of its radius. Dashed line I indicates the two-photon energy 2hu>i = 2.771 eV and solid line II 
corresponds to the scattered photon energy Tiuj s = 2hu>i — Tuoi t \ = 2.745 eV. The solid (dashed) arrows indicate the radii of the 
corresponding outgoing (incoming) resonances. Line III represents a Gaussian distribution over QD radii centered at R — 18 A 
and with FWHM equal to 40% (see text), (b) Frequency of the l p = 1 optical vibrons of a CdSe nanocrystallite as a function 
of its radius. 

FIG. 3. Plot of the electrostatic potentials associated to the first optical vibrons for different crystallite radii. Solid line: 
n p = 1; dashed line: n p — 2; dotted line: n p — 3; dot-dashed line: n p = 4. The equation u) npt i(R) = u>f with n v — 2, 3, and 4 
is fulfilled for R = 11.5 A, R = 16.2 A, and R = 21 A, respectively. We see that the electrostatic potential $n p ,i increases at 
the interface whenever co np ,i coincides with the Frohlich frequency u>f- 



FIG. 4. Hyper-Raman spectra of CdSe nanocrystallites embedded in glass for different radii and two photon energies: a) 
R = 16 A and 2Huji = 2.979 eV . b) R = 18 A and 2hoj l = 2.775 eV. c) R = 21 A and 27iw, = 2.559 eV. 

FIG. 5. Hyper-Raman intensity calculated for CdSe nanocrystallites with radius R = 18 A. The resonances are denoted by 
the labels IN or ON, where I (O) means incoming (outgoing) resonance with the exciton level N. Three approaches have 
been used: Solid line: Full matrix diagonalization including exciton effects. Dotted line: first order perturbation theory for the 
energy. Dashed line: free electron-hole model. 

FIG. 6. Hyper-Raman spectrum for 2huii — 2.771 eV of an ensemble of CdSe QD's with mean radius 21 A and a 40% size 
dispersion. The contributions of the QD's with different resonance radii are also shown. 

FIG. 7. Hyper-Raman spectra for an ensemble of QD's with a mean radius 18 A and a 40% size dispersion for different 
values of the laser energy. The inset illustrates the dependence of the maximum position and FWHM on the two photon energy 
2hcoi . The lines are a guide to the eyes. 

FIG. 8. Dependence of the HR lineshape with the mean radius R at 2hu>i — 2.771 eV. The size dispersion is 40% of the 
mean radius in all the spectra. 
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TABLE I. Values of the material parameters used for the 
numerical calculations. The same value P M (r„ ) has been 
assigned to all exciton (vibron) states. 



Parameter Value 

m e /m 0.12 a 

m h /m 0.45 a 

E g (eV) 1.865 b 

uj l (cm" 1 ) 213.1 b 

uj t (cm -1 ) 165. 2 C 

(3 L (lO^s" 1 ) 2.969 b 

(3 T (IO^is" 1 ) 0.002 

e a0 (CdSe) 9.53 d 

e aoc (CdSe) 5.72 c 

e bo o (Glass) 4.64 f 

r M (meV) 5 

r„ p (cm" 1 ) 2 

P 2 (eV) 20 g 



a Ref. H. 
b Ref. K 
c Ref. [M- 
d Ref. (32). 

c Calculated from the Lydanne-Sachs- Teller relation. 
f Ref. H. 

s Ref. f§. P 2 = 2|p c „| 2 /m 



TABLE II. List of resonance radii and corresponding reso- 
nant exciton levels (see text and Fig. 6). The incoming (out- 
going) resonances are labeled by IN (ON) 



R(A) 


18.0 


20.3 


26.0 


27.8 


28.0 


30.0 


(N,L) 


(1,0) 


(1,1) 


(2,1) 


(3,1) 


(3,0) 


(4,1) 
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